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We study the behavior of the defect and heat densities under sudden quenching near the quantum 
critical points in the two-dimensional Kitaev honeycomb model both in the thermodynamic and non- 
thermodynamic limits. We consider quenches starting from a quantum critical point into the gapped 
as well as the gap less phases. We choose points on the lines of anisotropic quantum critical points 
as well as different points of intersection of these lines as the initial points from where the quenching 
starts. We find that the defect and heat densities display the expected power-law scalings along with 
logarithmic corrections to scaling (or cusp singularities) in certain cases. In the vicinity of some of 
the intersection points the scaling behaviors change, indicating an effective dimensional reduction; 
the scaling behavior near these points depends on the number of critical lines crossed in the process 
of quenching. All the analytical predictions are also verified by numerical integration. 

PACS numbers: 64.70.Tg, 64.70.qj, 03.75.Lm 



I. INTRODUCTION 

A quantum phase transition (QPT) driven exclusively 
by quantum fluctuations at zero temperature is associ- 
ated with a dramatic change in the symmetry of the 
ground state of a many-body quantum HamiltoniarP^. 
It has been observed that quantum information theo- 
retic measures like the concurrence^, the entanglement 
entropy^KSl, fidelitjW^, etc., are able to capture the 
singularity associated with a QPT and they satisfy dis- 
tinct scaling relations close to the quantum critical point 
(QCP)(for review articles see |21|22j ). Also, the non- 
equilibrium dynamics of quantum critical systems and 
their connection to qua ntum information theory has been 
investigated extensivel y 23 * 24 * . The slow quenching dy- 
namics (defined by a rate) of a quantum system across a 
QCP and the scaling of the defect density and the heat 
density generated in the process have been major top- 
ics of research in this context! 25 " 27 !. The defect density, 
which represents the density of quasiparticles generated, 
and the heat density, or the excess energy above the new 
ground state, are expected to satisfy scaling relations 
given by the rate of driving and some of the exponents 
associated with the critical point across which the system 
is quenche d 27 1 28 1 (for review articles see |21|22|29j ) . 

In this paper, we are however interested in the scaling 
of the defect and heat densities for a sudden quench close 
to a QC P 18 * 19 * (see also [20]). In a sudden quench, when 
a parameter A in the Hamiltonian H of the system is 
changed abruptly, the wave function of the system does 
not have sufficient time to evolve. If the system is initially 
prepared in the ground state for the initial value of the 
driving parameter, it can not be in the ground state of 
the final Hamiltonian. Consequently, there are defects 
and excess energy in the final state. 

For a sudden quench of small amplitude from an initial 
value of a parameter A of the Hamiltonian to its final 
value A + 5, the defect density (n ex ) and the heat density 
(Q) can be calculated through the overlap between the 



ground state (^o(A)) of the initial Hamiltonian and the 
eigenstates (i/) n (\ + 5)) of the final Hamiltonian in the 
following way: 



l~ d I Wto(A)hMA + S))\ 2 (E n (X + 6) — E (X + 5)), 



where E n denotes the energy of the n'th excited state 
^ n (A + 5) of the final Hamiltonian. We note that n ex 
can also be expressed as n ex = (1/L d )(l — |(^o(A)|^o(A + 
(5))| 2 ), where \i/jq(\+5)) is ground state of the final Hamil- 
tonian. 

We consider a sudden quench of small amplitude A (i.e, 
S = A) starting from a QCP at A = 0; the defec t density 
can be related to the fidelity susceptibility 2 ^ Xf(A) at 
A as n ex = (l/L d )A 2 XF(A) where 



XF[ > ^ (E n (X) — E (X)) 2 ■ 



(2) 



Similarly, the heat density is related to the heat suscep- 
tibility xe as Q = (l/L d )A 2 X£;(A), where 



Xe 



(A) = ^l(^n(A)||f|^(A))| 2 



E n {\) - E (\) 



(3) 



which can be obtained by finding the overlaps 
(^o(A)|^ n (A + 5)) using adiabatic perturbation theory^. 
Both xf and xe exhibit interesting scaling behavior close 
to a QCP as discussed below. 

The scaling of n ex and Q follows from those of xf and 
Xe and is given by n ex ~ \ vd and Q ~ \ u ( d + z "> in the 
thermodynamic limit (L ^> X~ u ); here L is the linear 
dimension of a d-dimensional system and v and z are the 
correlation length exponent and the dynamical exponent 
associated with the corresponding QCP, respectively. In 
the opposite limit (L <C A _z/ ), the above scaling relations 
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get modified to n ex ~ \\\ 2 L 2 ^~ d and Q ~ \X\ 2 L 2 ^- d ~ z . 
It has also been predicted that the power-law scaling of 
n ex (or Q) is valid when vd (or v(d+z))< 2; otherwise the 
contribution coming from the low-energy modes becomes 
sub-leading and susceptibilities develop a cusp singularity 
at the QCflEl 

These scaling relations have been generalized to 
anisotropic quantum critical points (AQCPs)p21 (such as 
the ones appearing in the Kitaev model as discussed be- 
low), which have correlations length exponents v = z/m 
along m spatial directions and v = v± along the remain- 
ing (d - m) directions, and E n - E - A""*" = X u±z± , Q 
and n ex have the following power law scaling 



U ex ~ y\\rn+v ± {d-m) ^ 
Q )f\\m+i>±(d— m)+i/||«|| 



(4) 



for the thermodynamic (large system size) limit where 
A ^> Ln 1 ^ , L^ 1 ^ u± . Similarly, corresponding non- 
thermodynamic scaling relations can be written. 

We now discuss the motivation behind choosing the 
two-dimensional Kitaev modeP^ for the present work. 
To the best of our knowledge, this is the only integrable 
model in two-dimensions. Moreover, the phase diagram 
of the model (Fig. ^) has a gapless phase of finite width 
separated from the gapped regions by lines of AQCPeP. 
Although sudden quenching close to an AQCP has been 
studied earlier 30 , the Kitaev model provides us an op- 
portunity to investigate quenching into a gapless phase 
starting from an AQCP; interesting features arising due 
to the gapless nature of final state are emphasized in this 
work. Moreover, we study the quenching in the vicin- 
ity of the intersection points of the phase diagram where 
the model is effectively one dimensional. Interestingly, 
this effective "dimensional reduction" manifests itself in 
the scaling relations of the defect density even slightly 
away from these intersection points; one therefore ob- 
serves an interesting crossover behavior in the scaling as 
one approaches the upper intersection points along the 
boundary between the gapped and gapless phases. From 
the point of view of the predictions of adiabatic perturba- 
tion theory^, this model enables us to investigate power- 
law, logarithmic as well as cusp singularities in the heat 
susceptibility \e when the appropriate exponent is less 
than, equal to and greater than two, respectively. 

The paper is organized following way: in Sec. II, we 
discuss the model and propose the generic forms of the 
defect and heat density following a sudden quench. In 
Sec. Ill, we estimate the scaling relation for quenching 
into the gapped or gapless phase starting from an AQCP. 
In Sec. IV, the same is derived for the upper intersec- 
tion points while in Sec. VI quenching in the vicinity of 
the lower intersection point is investigated. In Sec. V, 
we present the numerical result that shows the crossover 
behavior of the scaling relation of the defect density dis- 
cussed above. Concluding remarks are presented in Sec. 
VII. 



II. KITAEV MODEL AND SUDDEN 
QUENCHING 

The two diemensional Kitaev model consisting of spin- 
1/2 's on a honeycomb lattice was initially proposed and 
solved by Kita ev in 200d^. The Hamiltonian for the 
mode l i jSTIMIMl 

j-\-l=even 

(5) 

where j and I respectively denote the column and row 
indices of a honeycomb lattice (see Fig. 1), and cr x -f' z are 
the spin- 1/2 operators (Pauli matrices) defined on the 
site j, /. We work with positive values of the J couplings. 
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FIG. 1: Schematic representation of the Kitaev model on a 
honeycomb lattice^ showing the bonds with couplings Ji, J2 
and J3. Mi and M2 are spanning vectors of the lattice. Sites 
'a*?' and '6^' represent the two inequivalent sites which make 
up a unit cell. The nearest-neighbor lattice spacing is taken 
to be unity. 

The hamiltonian can be diagonalized by the means of 
a Jordan- Wigner transformation^^ 

a 3,i = f n a ti ) a h f ° r even j + ^ 

\i= — 00 / 

a 'hi = ( II a ti) a li f or even j + 1, 



\t= — oo / 

b i,i = ( II °mW,J forodd i + 



6' 



\*= — OO / 

■3.1 - { n <jW for ° dd j+i > (6) 



which creates Majorana fermion operators a^, a'- ^ bjj 
and bj l that have their squares equal to 1, and anticom- 
mute with each other. 
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The sites j, I may alternatively be indexed by the vec- 
tors n = y/3im + (^i + |j)^2 which specify the centers 
of the vertical bonds of the honeycomb lattice, where i 
and j are the unit vectors along the horizontal and ver- 
tical directions respectively, and n\ and n 2 are integers. 
The Majorana fermions (a' n ) and b^ (b' n ) are located 
at the bottom and top lattice sites respectively of the 
bond labeled by n. The lattice vectors of the underlying 
triangular lattice are Mi = and M 2 ^+ |j- 

The nearest-neighbor lattice spacing is taken to be unity. 
In terms of the Majorana fermions, the Hamiltonian in 
Eq. (§ is given by 

n 

(7) 

where = i b f n a f H . The operators have eigenvalues 
±1, and commute with each other, and with H hence 
all the eigenstates of i^' can be labeled by specific values 
of Dfi. The ground state can be shown to correspond to 
Dfi — 1 for all n. 

The Fourier transforms of the Majorana fermions are 
given by 



-T 

k 



at e~ ik - f ' 

k 



(8) 



and similarly for a^, b^ and 6^. Here, N is the number 
of lattice sites, and the sum over k extends over half 
the Brillouin zone of the hexagonal lattice because of 
the Majorana nature of the fermiong^ 3 | 34 | 36 i The full 
Brillouin zone is taken to be a rhombus with vertices 
lying at (k x ,k y ) = (±2tt/>/3,0) and (0,±2tt/3); half the 
Brillouin zone is given by an equilateral triangle with 
vertices at (k x ,k y ) = (27r/-\/3, 0) and (0,±27r/3). 

For = 1, the Hamiltonian can be diagonalized into 
the form 



(9) 



where Hj: can be written in terms of Pauli matrices as 



where 
and f3r 



a k al + h a ^ 

2[Ji sin(£ • Mi) - J 2 sin(k • M 2 )], 
2[J3 + J\ cos(k ■ Mi) + J 2 cos(k • M 2 )]. 

(10) 



The energy spectrum of consists of two bands with ener- 
gies given by^ 



Et = 



± 



where 
/% = 2[J 3 



Ji cos( — k x 



2 ky) 



a n = 2[J 1 sm( — k x - -k y ) 



J 2 cos( — k x 



(11) 



J 2 sin( — fca; 




(b) 
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FIG. 2: (color online) (a) Phase diagram of the Kitaev model 
in 3-D. The interior of the green wedge is the gap less phase, 
where the couplings satisfy the triangle inequalities given by 
Ji < J 2 + ^3, J 2 < J3 + J\ and J3 < Ji + J2. The yellow 
plane is Ji + J2 + J3 = 2, the projection onto which gives (b). 
The green wedge of the gapless phase extends indefinitely but 
has been truncated for clear representation in the figure, (b) 
Projection of the phase diagram onto the plane J\ + J2 + J3 = 
2. The gapless phase is the region inside the inner equilateral 
triangle. All the points on the boundaries between the gapless 
phase and the gapped phases are AQCPs (see Appendix A). 



The energy gap Et 



Ej* vanishes for specific values of 



k when J3 < J\ + J 2 and J\ < J 2 + J3 and J 2 < J3 + J\ 
giving rise to a gapless phase of the model. The gapless 
and gapped phases of the model are shown in Fig. [2] The 
phase boundaries between the gapless and gapped phases 
are given by J\ + J 2 = J3, J 2 + J3 = J\ and J 3 + J\ = J 2 . 
On these boundaries, away from their intersection lines, 
we have anisotropic quantum critical points, i.e., the dis- 
persion varies with linearly with k along one direction 
and quadratically along another distinct direction (see 
Appendix A). 

The ground state is given 



1*0) 



n 



I (4 



|*>. (13) 



where the product runs over half the Brillouin zone, 
a/e, a f k , b' k are the Fourier transforms of Majorana 
fermion operators used, and 



-I 



(14) 



Excited states are produced by exciting modes corre- 
sponding to k' \—k' to the upper band and are given by 



n 



1 



q {a ^ +e ^, h ^ ){a n +ih n i)) 



U4-e^bl)(a'Uib'h 



|$). (15) 



ky)}. 



(12) 



We consider the coupling J3 as the quenching parame- 
ter. It's initial and final values are given by J3 = J^ c + A 
and J3 = Jsc + A + S respectively, where Js c is the value 



4 



of J3 at the QCP and 5 is the small quench amplitude. 
The defect and heat densities can be calculated from the 
overlap of states as given in Eq. (§ can be expressed as^ 



7 max 



^ Jtt/L h 

± 5 > r x r~ 

^ L/l Lil %(A)2|£7g(A + <y)| 



2 dkydlvgr j 



k 



dk y dk Xl (16) 



where Er is given in Eq. (11). Here the integrals run 



across half the Brillouin zone. The wave vector k is mea- 
sured with respect to the wave vector for minimum en- 
ergy. (Note that the 1/L 2 term in Eq. Q gets cancelled). 
For a quench starting at the QCP (A = 0) and ending at 
some small value A (i.e., 5 = A), the integrand can be 
expanded about A, giving the values of n ex and Q to 
leading order in A as 



■1 max imax O 



tt/L Jtt/L 

„max imax 



tt/L Jtt/L 



dky dk x 



;dk y dk x (17) 



In the remainder of this work, we shall use expansions of 
and /3g about the minimum energy points to evaluate 
the integrals for analytic calculations, and directly eval- 
uate the whole integrals over half the Brillouin zone for 
numerical calculations. 



A. Defect density 

We first consider quenches starting in the AQCP A 
and ending in the gapped phase (J3 = 1 + A) with A <C 1. 
The analysis presented in this section also applies to the 
other AQCPs under quenching of J3 (appendix A). In 
this limit, one can neglect the k 2 term in the expansion 

and k y to 



of j3 in Eq. (18) and extending the limits of k x 



infinity one gets 



9A 2 

7T 2 



J TT I L J TT 



k y 



dky dkgg , 



(20) 



where E± = yj(9k% + (§/c 2 ± 2A) 2 ). Using the results 

presented in reference [3Q for the scaling of the fidelity, 
one can readily show that the defect density scales as 
A 3 / 2 in the thermodynamic limit and n ex ~ A 2 L X / 2 in the 
other limit which is consistent with the expected scaling 
of the defect density at an AQCP (as given in [l| with 
d = 2,m = 1,z/_l = 1 and ^ = I/2P. 

For quenches starting from the AQCP into the gapless 
phase, we use Eq. (19) and employ the expansion about 



the gapless points retaining only the leading order in k x 
and k v 



9A 2 r 

^ 2 L 



/ L Jtt/L 



-dk y dk x . (21) 



One should note here that a quench of amplitude A ex- 
cites modes up to k\\ = k x ~ A^n = y/X and k± = k y ~ 
X u± = A; we then propose a scaling form 



III. SUDDEN QUENCHING INVOLVING A 
SINGLE CRITICAL LINE 



1 



oc 



-/ 



1 



1 



L\\L± \ Li 1 -s/X' L±X 



(22) 



We first consider the point J\ = J 2 = 1/2, J3 = J\ + 
J 2 = 1 (point A in Fig. (|2|). The energy gap is zero at 
the four corners of the Brillouin zone. Now, if J3 is offset 
slightly into the gapped phase, it leads to small energy 
gaps at these points. Expanding and /?£ in terms of 
the deviations k x , k y from the zero energy points, we get, 
for J 3 = 1 + A 



3ky , 

3 2 9 2 

4^+4% 



2A. 



(18) 



where / is the scaling function. We assume the thermo- 
dynamic limit L\\y/X ^> 1, L^A ^> 1 and L\\ = L± = L; 



the first argument of the scaling function in Eq. (22) is 



negligible in comparison to the second argument. We 
then have 

n ex oc A^/(L_lA), (23) 

where /(Lj_A) = /(0, l/Lj_A). To find the form of /, we 
evaluate the relevant integral for n ex (see Appendix B), 
with the appropriate cut-offs of y/X and A for k x and k y 
respectively and Lm = L± = L, which gives 



On the other hand, if J3 is changed by a small amount 
into the gapless phase, the gapless point of the dispersion 
is shifted from the corners. In the gapless phase the ex- 
pansion can also be done about the new shifted gapless 
points, which is more convenient for certain calculations 



h = \ki - 



^ - V6Xk x . 



(19) 



A 3 / 2 ln(A£). 



(24) 



Eq. ( 24 ) presents the exact analytical result for the de- 



fect density for quenching into the gapless phase and 
shows that the logarithm that appears in the expression 
for tIq X is a function of a dimensionless combination of A 
and L (more precisely L±). Numerically, n ex ~ A 3 / 2 In A 
(with L fixed) and n ex ~ InL (for fixed A) (see Fig.Q); 
we attribute these logarithmic corrections to the scaling 
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A 3/2 n ex 






* 3/2 n ex 


In L 






In A 



FIG. 3: Variation of \~ 3 / 2 n ex with In A in the thermodynamic 
limit for a quench into the gapless phase obtained numeri- 
cally for Jx — J 2 — 1/2, J3 = 1, L = 10000; n ex varies as 
A 3/2 ln(AL). (Inset) Variation of A" 3/2 n ex with InL in the 
thermodynamic limit for a quench into the gapless phase ob- 
tained numerically for Ji = J2 = 1/2, J3 = 1, A = 0.002; n ex 
varies as A 3 ^ 2 ln(AL). 



to the gapless nature of the phase into which the system 
is quenched to. 

We note that in the nonthermodynamic limit (L « 
A -1 ), the integrand does not have any singularities and 
one gets identical scaling as that for quenching into the 
gapped phase. 









A 



FIG. 4: Variation of xe = Q/A 2 with A in the thermodynamic 
limit obtained numerically for Ji = J 2 = 1/2, J 3 = 1, L = 
10000. It can be seen that this quantity does not diverge or 
go to zero at A = 0; Q varies as A 2 . There is a cusp singularity 
in xe at A = 0. 



in Eq. (18) (as was done while calculating n ex ), thereby 
focussing only on the low energy modes, one gets the re- 
lation xe °° l^l 1 ^ 2 ? this scaling gets modified because of 
the contributions of the high energy modes. It is straight- 
forward to show that the above interesting cusp-singular 
behavior can be seen at any other AQCP of the phase 
diagram. 



IV. SUDDEN QUENCHING INVOLVING THE 
UPPER INTERSECTION POINTS 



B. Heat density 

To calculate the heat density following a small sud- 
den quench of amplitude A into the gapped (or gap- 
less) phase starting from the point A (see Fig.([2|), we 
use Eq. (17) and employ the coordinate transformation 



k = ^(3/4)/c 2 + (9/4) & 2 , k y = k y . For the heat suscep- 
tibility xe, we have 



Xe 



T2 f°° 
V3tt 2 Jo Jo 



kyk dkydk 



(9^ + (p + 2A) 2 )ly / A 



9p 

(25) 

This gives xe ~ 16/(7p\/3) at A = 0, indicating that there 
is no divergence at A = even in the limit of infinite 
system size. The leading order term in the scaling of 
Q = \ 2 Xe, is thus A 2 instead of A 5 / 2 (which is sub- 
leading for A > 0) as expected from the scaling relations 
discussed in Sec. Uf^. 

The most interesting feature associated with xe is that 
it has a cusp singularity at A = as shown Fig.Q. The 
heat susceptibility xe nearly equals 16/ (^V^) for A < 
(i.e., in the gapless phase) and 16/(^^/3) — 72y/\/ (7r 2 \/3) 
for A > (gapped phase) forming a cusp at A = 0. 

The origin of this cusp singular behavior of xe stems 
from the following exponent relations; in the present case, 
v\\m + v_\_(d — m) + v\\z\\ exceeds two so that any power- 
law scaling predicted from adiabatic perturbation the- 



ory is sub- leading 18 . If one ignores the ky term in (3 



Here we consider the intersection point J2 = J3 = 
1, J\ =0 (point B of the phase diagram Qb)); the anal- 
ysis for the other intersection point J\ = J3 = 1, J2 = 
is completely equivalent); at this point the coupling J\ 
vanishes, effectively giving de-linked spin chains run- 
ning along the M2 direction. Expanding a& and /3k 
around the wave vector for the minimum energy with 
J3 = 1 + A, we find = 2k and /3g = k 2 + 2A, where 
k = (y/3/2)k x -\-(3/2)k y . Therefore there is always a gap 
in the spectrum when A ^ 0. Moreover, one effectively 
arrives at a one-dimensional dispersion with v = z = 1. 



A. Defect density 

To calculate the defect density for a sudden quench of 
amplitude A starting from the point B, we use Eq. (17) 
to obtain 



n ex -^\' 



I 

J 7T 



/L (4fc 2 + (fc 2 + 2A) 2 ) : 



-dk. 



(26) 



Using the thermodynamic limit of A 3> 1/L and rescaling 
k = \k', one finds 



1 



A 



k r 



4tt 2 J (k' 2 + l) : 



-dk!. 



(27) 



Thus n ex ~ A, which is consistent with the scaling rela- 
tion of n ex when d = 1 and v = 1. This is numerically 
confirmed in Fig.([5|. 
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FIG. 5: Variation of n ex with A in the thermodynamic limit 
obtained numerically for J2 = J3 = 1, Ji = 0, L = 10000; 
Hex varies as A. (Inset) Variation of n ex with L in the non- 
thermodynamic limit obtained numerically for J 2 — J3 — 
1, Ji =0, A = 0.0001; n ex varies as L. 



In the non- thermodynamic limit of A <C 1/L, the A in 
j3 can be ignored to give 



;X 2 



k 2 



" 7T 2 J„ /L (4k*+k*)* 

Thus, to leading order, 

1 



dk. 



nex — A L ~ A L. 



47T 3 ' 



(28) 



(29) 



This result is also expected from scaling for a one- 
dimensional QCP with v — 1. 



FIG. 6: Variation of xe — Q/A 2 with In A in the thermody- 
namic limit obtained numerically for J2 = J3 = 1, J\ — 0, 
L = 10000; xe varies as ln(l/A) = — In A. (Inset) Variation 
of Xe — Q/A 2 with InL in the non-thermodynamic limit ob- 
tained numerically for J 2 = J3 = 1, Ji = 0, A = 0.0001; xe 
varies as InL. 



anisotropic critical line in Fig.fj^b)), this value grows and 
diverges logarithmically as the horizontal distance from 
the point B vanishes. 



SUDDEN QUENCHING INVOLVING TWO 
CRITICAL LINES 



B. Heat density 

To calculate xe in the limit of infinite system size we 
retain only lowest order terms in (17) so that 



Xe 



16 

^2 



Jo 



(4k 2 + 4A 2 ) 3 / : 



-dk. 



(30) 



with k max ^> 1/L. Hence, to leading order, \e — 
2 / 'tt 2 (ln(2k max I A)) so that xe ~ In A -1 and Q ~ 
A 2 In A -1 . The heat susceptibility thus diverges logarith- 
mically at these points. 

For the non-thermodynamic limit, 



Xe 



16 

^2 



k 2 



ir/L 



(4k 2 + k A f/- 



-dk, 



(31) 



(2/7r 2 )ln(L/7r). Thus 



giving, to leading order, \e 
Xe ~ InL and Q ~ A 2 InL. 

We thus find a logarithmic divergence of xe which is 
confirmed numerically in Fig.(|6|. It is to be noted that 
in this case (d + z)v = 2; for (d + z)v < 2 (> 2) one 
has a power-law (cusp) singularity as already seen in the 
previous sections^. This case happens to be the marginal 
case, where one encounters logarithmic singularities. We 
have already shown that xe has a finite value at the point 
A; as one approaches the intersection point B along the 



We now consider a sudden quench from J3 = 1 + A to 
J3 = 1 — A in the vicinity of the intersection point B; we 
set J\ = e, J2 = 1 — e, J3 = 1 with e << 1. As can be 
seen from the phase diagram Qb)), such a quench will 
make the system cross one or two critical lines depending 
on whether A < 2e or A > 2e, respectively. Numerical 
investigation of the scaling of the defect density (with A) 
in the thermodynamic limit shows that: (i) when A — >• 
0, e — ^ 0, A < 2e, the defect density scales as A 3 / 2 and (ii) 
when A — » 0, e — >• 0, A > 2e, the defect density scales as A 
(see Fig. (ffl 



The above numerical result can be justified as follows: 
when A < 2e, the quench is between a gapped and gapless 
phase separated by an AQCP, and the A 3 / 2 scaling is 
expected from the scaling relation. The crossover to the 
scaling n ex ~ A occurs because the system crosses two 
anisotropic critical lines (starting from one gapped phase 
and reaching the other) and and the scaling eventually 
must approach that of a quench through the intersection 
point B as e — >• 0. 

We therefore find an effective "dimensional reduction" 
close enough to the intersection point so that the system 
crosses two anisotropic critical lines. It should be noted 
that a similar dimensional reduction is seen in the scaling 
of the defect density in the final state of the Kitaev model 
following a slow quencfP^. 



7 



ln(n ex /A 2 ) 


A=2s 




XlnA 



FIG. 7: Variation of n ex /A 2 with In A in the thermodynamic 
limit obtained numerically for J\ — 0.002, J 2 = 0.998, J3 = 1, 
L = 10000; n ex /\ 2 varies as A~ 1/2 for A < 2e = 0.004 and 
A -1 for A > 2e = 0.004. There is a clear change of slope at 
A = 2e. 



VI. SUDDEN QUENCHING INVOLVING THE 
LOWER INTERSECTION POINT 

We consider sudden quenches starting at the lower in- 
tersection point C (Ji = J 2 — 1, J3 = 0) into the gap- 
less phase above it. Following similar calculations, as for 
the point B, one finds in the non-thermodynamic limit 
n ex ~ A 2 L and Q ~ A 2 InL; this are identical to the 
scaling relations obtained for quenching from point B. 

In the thermodynamic limit, on the other hand, a sim- 
ilar analysis as given in Sec. Ill A for the quenching into 
gapless phase leads to the scaling n ex ~ Aln(AL); it 
should be noted that the scaling for tiqx is a that of the 
intersection point B along with a logarithmic correction 
arising from the gapless phase. Numerically one finds 
n ex ~ A In A and n ex ~ In L (see Fig. J8l) 



the quenching scheme in which only J3 is varied. 



VII. CONCLUSION 

We have studied the behavior of defect and heat den- 
sity under sudden quenches of the parameter J3 in the 
Kitaev model. The main results of the paper are the 
following: for quenching from an AQCP to the gapless 
phase, we have found an exact scaling form for n ex and 
Q. The susceptibility \e shows a cusp singularity at the 
AQCP confirming the prediction of adiabatic perturba- 
tion theory. Moreover, n ex has a logarithmic correction 
which carries the signature of a gapless phase. For a sud- 
den quench starting from the upper intersection point, 
we retrieve the scaling relations of the equivalent one- 
dimensional system for the defect density while there is 
an additional logarithmic correction in xe- What is more 
interesting that one observes a crossover in the scaling of 
the defect density depending on number of critical lines 
crossed in the process of quenching; very close to this 
intersection point, the system crosses two critical lines 
even for a quench of small amplitude and hence one finds 
a scaling expected for a one-dimensional system. For 
quenching into the gapless phase starting from the lower 
intersection point, one finds an additional logarithmic 
correction to the defect density only in the thermody- 
namic limit. All our analytical studies are supported by 
numerical calculations. 
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India via the KVPY fellowship. AD acknowledges CSIR, 
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FIG. 8: Variation of X~ 1 n ex with In A in the thermodynamic 
limit for a quench into the gapless phase obtained numerically 
for J\ — J2 — 1, J3 = 0, L = 10000; n ex varies as A -1 ln(AL). 
(Inset) Variation of X~ 1 n ex with InL in the thermodynamic 
limit for a quench into the gapless phase obtained numerically 
for Ji = J 2 = 1, J3 = 0, A = 0.002; n ex varies as A -1 ln(AL). 



Appendix A: Anistortopic Quantum Critical Points 

The phase boundaries between the gapless and gapped 
phases are given by J\ + J 2 = J3, J 2 + J3 = J\ and 

+ = ^2- On these boundaries, away from their 
three intersection points, we have anisotropic quantum 
critical points, i.e., the dispersion varies with differ- 
ent powers of k along different directions. Away from 
the intersection points of the lines, for J\ + J2 = J3, 
the dispersion goes to zero at distinct points given by 
(k x ,k y ) = (±2tt/\/3,0) and (k x ,k y ) = (0,±2tt/3). This 
happens at (k x ,k y ) = (±7r/\/3, T 71 "/^) for J2 + J3 = J\ 
and (k Xj k y ) = (±ir/y/3, ±7r/3) for J3 + J\ = J2. If 
and /?£ are expanded in terms of the deviations k x , k y 
from these points, then in general, 



In the thermodynamic limit, the heat density is found 
to scale as Q ~ A 2 ln(l/A), which is the same scaling as 
that for the point B. There is no additional log correc- 
tion from the gapless phase to the scaling of the heat 
density. It is to be noted that the point C is not equiva- 
lent to the point B because of the asymmetric nature of 



a\k x 

m2- 



- a 2 k y = fci, 
fak% + b 3 k x k y , 



(Al) 



where 61, 62 > 0. varies linearly and /3g varies quadrat- 
ically along any direction, due to the sine and cosine 
terms in a and j3 respectively. There is no constant term 
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in the quadratic polynomial for /3 k as the energy gap must 
vanish at k = 0. We now define k 2 = k x cos (j) + k y sin <\> so 
that there is no k\k 2 term when /3k is expressed in terms 
of k\ and k 2 . This gives 



tan ( 



2aib 2 - a 2 b 3 
ai6 3 - 2a 2 6i ' 



(A2) 



which always has a solution in 



-§,§]. We then have 



(A3) 



where 



c 2 



bi sin 2 (f) -\- b 2 cos 2 (f) — b3 sin cos ( 



(ai sine 



a 2 cos ( 



&ia 2 + 62a 2 — b 3 a\a 2 



(ai sin ( 



(22 cos ( 



(A4) 



Thus, the dispersion a 2 + /3 2 will vary linearly along 
fci when &2 = and quadratically along &2 when fci = 0. 
Hence we have an AQCP with v\ = 1, v 2 = 1/2. If we 
offset J3 by an amount A, 



/3g = cik 2 



c 2 k\ 



2A. 



(A5) 



For the first line c 2 > and hence the gap vanishes for 
A < 0. For the second and third lines, c 2 < and the 
gap vanishes for A > 0. These can be seen from the phase 
diagram, where the gapless phase is below the first line 
and above the second and third lines. The gap in the 
spectrum at k\ = k 2 = goes as A and hence v\Z\ = 
^2^2 = vz = 1. Since k\ and k 2 are linear functions of 
k x and k y , the Jacobian in the transformation of the 
coordinates is simply a constant. Hence, no extra factors 
that are functions of k appear in any integral over k space 
when changing coordinates. 

The expansions of a% and /3% in terms of k\ and k 2 
give, using Eq.(ll) and Eq.(17), 



A 



k\dk\dk 2 



Q 



L/lL/l (fc? + ( Cl fc? + c 2 fc 2 2 ±2A)2)2' 

2 f°° f 00 k\dk^dk 2 

L/lL/L (fc? + ( Cl fc?+C 2 fc 2 2 ±2A)2)3/2' 



which determine identical scaling behavior of n ex and Q 
at all AQCPs. 



Appendix B: Calculation of Defect Density for a 
Quench into the Gapless phase 



We have, for the defect density, 
9A 2 f x kl 



tt/L Jtt/L 



(9fc2 + 6Afc2) 



2 dk y dk x . (Bl) 



Integrating with respect to k y , we obtain 



9A 



2 r VX 



f (G 1 (X)-G 1 (y) + G 2 (X)-G 2 (y))dk x 

Jtt/L L L 

(B2) 



where 



G 1 (ky) = 
G 2 (ky) = 



6(3*2 + 2k%\) 



6V3k x V2X 



(B3) 



Integrating with respect to k x gives 



q \2 » 



7T^ 



(B4) 



3=1 



where 



T 2 



tan x (y |) ^ -0.685 

6vWa ~ 6vWa 



6vWA 6V6LA 

tan -i(4^) 



T 7 



6>/6\/A 12vWA 
tan-^y^fv^) 1 
18 

1.085 



12\/6>/A 

7T 



12LM 



7rln(y / §v / A) 
12VWA 



(B5) 



The most significant terms are and Tg because of the 
large In terms. Hence n ex ~ (3/(4ttv / 6))A 3 / 2 ln(LA/7r). 
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